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ABSTRACT 

This paper formulates time-dependent, force- free, degenerate electrodynamics as a 
hyperbolic system of conservations laws. It is shown that this system has four char- 
acteristic modes, a pair of fast waves propagating with the speed of light and a pair 
of Alfven waves. All these modes are linearly degenerate. The results of this analytic 
study can be used in developing upwind numerical schemes for the electrodynamics of 
black hole and pulsar magnetospheres. As an example, this paper describes a simple 
one-dimensional numerical scheme based on linear and exact Riemann solvers. 

Key words: black hole physics - pulsars:general - magnetic fields - methods: numerical 



1 INTRODUCTION 

In magnetospheres of pulsars and black holes the electromagnetic field is so strong that the inertia and pressure of the 
plasma can be ignored. As a result the Lorentz force, F a/ 3 J a , vanishes and the transport of energy and momentum is entirely 
electromagnetic ( Gtoldreich fe Julian 196^ ; Blandford & Znajek 1977). This justifies the name "force-free" to describe the 



electrodynamics of pulsars and black holes. However, the electrodynamics of the magnetospheres is rather different from 
that in a vacuum which is obviously also force-free. Indeed, the magnetospheric plasma is plentiful enough to support strong 
electric currents and screen the electric field ( Goldreich fc Julian 1969} Blandford & Znajek 1977). An electromagnetic field 
satisfying this condition is called "degenerate" (Macdonald & Thorne 1982). Thus, the electrodynamics of pulsar and black 
hole magnetospheres is force-free, degenerate electrodynamics (FFDE). 

So far, theoretical studies of pulsar and black hole magnetospheres have assumed a steady-state, indeed the properties of 



FFDE as a system of time-dependent equation have not been studied systematically. It is only recently that Uchida (1997 



has developed a rather elegant theory of this system in which the electromagnetic field is described in terms of two scalar 
functions, called "Euler potentials". However, this formulation is not particularly suitable for numerical analysis because 
its basic equations, when written in components, involve mixed space and time second order derivatives. In this paper, we 
present another formulation in which the evolution equations have the form of hyperbolic conservation laws. We also describe 
a simple one-dimensional, upwind numerical scheme based on these analytic results. This scheme can easily be generalized to 
mu ltidimensions and curved space-time. In fact, we have already c onstructed a 2D scheme f or the Kerr space-time and used 
it ( Komissarov 2001 ) to study the Blandford-Znajek mechanism ( Blandford fc Znajek 1977 ) for the extraction of rotational 
energy from black holes. 



2 EVOLUTION EQUATIONS OF DEGENERATE ELECTRODYNAMICS 



Since far away from the relativistic "star" the inertia of particles is bound to become important, the approximation of ideal 



relativistic magnetohydrodynamics (RMHD) could be considered as more satisfactory in general (Phinneyl982). This suggests 
one way of deriving the time-dependent equations of FFDE, namely by considering the low inertia limit of RMDH. We shell 
adopt this approach because it reveals the close connection between these two systems. However, the time-dependent equations 
of FFDE can also be derived entirely within the frame of electrodynamics in a rather obvious way. Besides that, some of the 
microphysical conditions of MHD, like the small mean free path of plasma particles, do not have to be satisfied in FFDE. 



In the covariant form the equations of RMHD are as follows ( Lichnerowicz 1967; Anile 1989) 
Maxwell's equation: 
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V Q *F a0 = 0, (1) 
Energy-momentum equation: 

V a T aP = 0, (2) 
Continuity equations 

V Q (n w w a )=0, (3) 

where T al3 is the total stress-energy-momentum tensor, *F al3 is the dual tensor of electromagnetic field, nu) is the proper 
number density for species i and u a is the fluid 4-velocity. The total stress-energy-momentum tensor 

rpOL/3 _ rpaf3 . rpa/3 f,\ 

is the sum of the stress-energy-momentum tensor of matter, T?^-. , and the stress-energy-momentum tensor of electromagnetic 
field 

T$ = F^F 9 " - \{F^ v F^)g af \ (5) 

where g a/3 is the metric tensor and F a /3 is the electromagnetic field tensor. In this paper we assume that Greek indices run 
from to 3 and Latin indices from 1 to 3. 

These differential equations (^-[]) are supplemented by equations of state and the perfect conductivity condition 

F^u" = 0. (6) 

This condition makes the second group of Maxwell equations redundant. From ^ one obtains 

F -*F = 0, (7) 

F ■ F > 0, (8) 

where F ■ F means F llv F^ 1 ' . Just like classical MHD, RMHD is a hyperbolic system of conservation laws and has waves of 
four types - entropy, fast, Alfven, and slow (e.g Anile 1989; Komissarov 1999). 

If the contribution of matter to the total tensor of stress-energy momentum is very small then a perturbation technique 
suggests itself. The zero order equations constitute the system of FFDE: 

V M •f" = 0, (9) 
V„T$ = 0, (10) 
From (^|^|) one has 

v„r$ = -f^j\ (ii) 

where 

J" = V^F"" 

and, thus, (^) requires the vanishing of the Lorentz force. Since equations (PJlcl) do not involve u a , the perfect conductivity 
condition should be replaced with conditions (QJ^), which express the degeneracy of electromagnetic field in FFDE. They 
ensure the existence of observers who detect only magnetic field. Let a 11 to be the velocity of such observer. It easy to see that 

F^a v = 

and, thus, a" is a zero eigenvector of F The second zero eigenvector can then be introduced via 
In terms of a M and b M 

F^v = e fl „ Q/3 a a b /3 , (12) 
and 

*F> iv =26 [ V'. (13) 

From the last equation it follows that the unit space-like vector c M orthogonal to a M and 6 M is a zero eigenvector of *F M1/ . 
The second zero eigenvector, also space-like and orthogonal to a M and , is given by 
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In terms of these vectors 

*F M „ = e^ a pc a d 13 , (14) 
and 

F M " = 2c 1m cT 1 . (15) 

From (|t]) it follows that there are only five independent components of -F M „. In components, (0) splits into three evolution 
equations and one differential constraint on the initial data (XTjB 1 = 0). Thus, if the system (Q)-([lo]) is self-consistent then ( |lo| ) 
has only two independent components. Indeed, equations (|ll],[l2|) show that that in the limit of degenerate electrodynamics, 
the vector V M T yM always belongs to the two-dimensional vector space generated by and d M . 

The evolution equations of the system ([],|l^) have the form of conservation laws. In the following section we shall show 
that this system is also hyperbolic. 



3 HYPERBOLICITY OF DEGENERATE ELECTRODYNAMICS 

In order to demonstrate the hyperbolicity of FFDE, it is sufficient to study its one-dimensional equations in a locally pseudo- 
Cartesian system of coordinates. In such a frame, the electric and magnetic fields are defined via 



E i = F i0 , B*=^ K F jh , (16) 
where e 1 ^ is the 3-dimensional Levi-Civita symbol. Then 

T H = _( E i E 3 + B l B 3 ) + ^-{E 2 + B 2 ), 
and conditions ([?],[]) read 

E ■ B =0, (17) 

B 2 -E 2 >0. (18) 

Notice that we use units such that the speed of light and 7r do not appear in the equations. 

Like in MHD, the number of evolution equations and dependent variables for the one-dimensional system is reduced by 
one. Indeed, if d/dx 2 = d/dx 3 — then x°- (t-) and the a; 1 -components of (|^) require 

B\ — const. 

Thus, the original system of FFDE is reduced to four independent evolution equations for four dependent variables, B2, S3, 
and two components of the electric field. The selection of these two components depends on the direction of magnetic field. 
For example, if B\ — then one of them must be E\. As a result, we have to study three different sets of ID equations to cover 
all possible cases. One can avoid such complications by relaxing the degeneracy condition (|l7]) and modifying the equations 
in a way such that this condition is satisfied automatically if it is satisfied by the initial data. For example, one can select all 
three space components of the energy-momentum equation and add the term Bid(BiE l )/dt to obtain 

d_F^ dF* d(B^) 

m + dxl +B t gt - U. (19) 

Together with the two remaining components of (|^) these equations constitute what we shell call the "augmented one- 
dimensional system of FFDE". It has five independent equations for five components of F vli . Obviously, solutions of this 
augmented system satisfying the constraint ( |l7| ) are also solutions of the original system of FFDE. Moreover, contraction of 
@ with B { gives 

ox L at 

This result shows that if B 1, E i = at t = then all time derivatives of B l Ei at t = also vanish and, therefore, B 1 E i = at 
t > as well. 

All characteristic waves of the FFDE must be present in the augmented system because of the way it has been constructed. 
It is also expected to have one additional nonphysical wave due to higher number of evolution equations. This extra wave can 
easily be identified because across this wave EiB 1 does not have to stay constant. 

In vector form the augmented system reads 
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A dU „dU „ 

where 

U={B 2 ,B3,Ei,E2,E 3 ) t , 

is the vector of dependent variables, 



^4 = 
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(20) 



(21) 



(22) 



c 



Given ( |rij , |ilj[ ) the eigenvalue problem 

(A - /iC) • r = 

has the following solutions 

nf = ±1, 



(23) 



± S 3 S 2 - S2S3 ± JB\{B?- -S 2 ) 

Ma = 



i? 2 



/in = 0. 



(24) 

(25) 

(26) 
(27) 



It is easy to verify that in the limit of FFDE the wave speeds of the fast waves of relativistic MHD ( Anile 198£ ; Komissarov 
19980 reduce to ^ and the wave speeds of Alfven waves of relativistic MHD reduce to . Hence, fj,„ corresponds to the 
nonphysical wave of the augmented system. 

If Si 7^ then the eigenvectors of the fast modes are 

rf = (-T)f, </, 0, nfCf, nfvf) , 
where 



ri f = E s + iJ,fB2, Cf /•:■ /<./>':. 

It is shown in Sec. 4 that these waves have the same properties as linearly polarised electromagnetic waves in vacuum. 
The eigenvectors of the Alfven modes are 

(vl + C|) ± _ ± t 

1 



where 

Tja = S3 + fttB 2 , <a = E 2 



Ma S3, 



(28) 
(29) 

(30) 
(31) 



In some applications it might be useful to determine the left eigenvectors as well as the right eigenvectors. Unfortunately, 
the author was unable to find a concise form for the left eigenvectors and they are therefore not presented here. 

Like MHD, the system of FFDE is not strictly hyperbolic, that is it allows multiple eigenvalues under certain conditions. 
Such degenerate^ cases require separate treatment because the eigenvectors given by (^8|) and (^o|) become either singular or 
linearly dependent. 

(i) Provided 

E t = -i 1 xB, (32) 
where i± is the unit vector along the £i-axis and Et is the tangential component of electric field, one has 



* In this paper three different kinds of degeneracy are encountered: 1) the degeneracy of electromagnetic field itself (see Sec.l), 2) the 
multiplicity of eigenvalues of the Jacobean matrix (Sec. 3), and 3) the so-called linear degeneracy of hyperbolic waves (see Sec. 4). 
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Ma = M/ = +1- 

and the corresponding eigenvectors form the two-dimensional vector space generated by 

(0, 1, 0, 1, 0)* and (-1, 0, 0, 0, 1)'. (33) 

Obviously, vacuum electromagnetic waves propagating in the positive direction have the same eigenspace. Thus, in this limit 
the right Alfven wave also behaves like a linearly polarised electromagnetic wave in vacuum. Similarly, if 

E t = nxB, (34) 

one has 

Ma = M/ = "I, 

and the corresponding eigenvectors form the two-dimensional vector space generated by 

(0, 1, 0, -I, 0)' and (1, 0, 0, 0, 1)'. (35) 

This eigenspace is the same as that of vacuum electromagnetic waves propagating in the negative direction. We notice that 
conditions (^) and ( |3^ ) cannot be satisfied simultaneously unless E = and Bt = 0. 
(ii) If Si = then 

Ma = Ma = Ma = "gg = ~ "g" (36) 

with a two-dimensional vector space of eigenvectors. One can use 

r = (l, 0, ~ jB2( g°~ 1) , 0, (37) 
and 

r = (0, 1, ~ j33( ^~ 1) , ^, 0) 4 (38) 
as the basis of the eigenspace if E\ 7^ or 

r = (0, 0, 1, 0, O) 4 (39) 
and 

r = (£3, -Ba, 0, -ti a B 2 , -ytaBzf (40) 
if E 1 = 0. 



(41) 



4 PROPERTIES OF WAVES OF DEGENERATE ELECTRODYNAMICS 
4.1 Fast waves 

From ( pg| ) one obtains the following system of differential equations describing the fast simple wave: 

dB 2 _ dB 3 _ dE 1 _ dE 2 _ dE 3 

-(E 3 + n f B 2 ) ~ E 2 -HfBz ~ ~0~ ~ fi f (E 2 - fi f B 3 ) ~ ^ f {E 3 + ^ s B 2 y 

where fj,f is either jiti" = +1 or fxj — —1. These are easily integrated to obtain 

Ex = const. (42) 

Tjf — E 3 + )ifB 2 — const, (43) 

C/ =Ea- VfB 3 = const, (44) 

E i B i = const, (45) 
If we introduce the vector 

*> = (o,Cf,vfY 

then 

dE t \\t f and dB t ±t f , (46) 

where Bt and Et are the tangential components of the fields. Thus, these are transverse waves and they have the same 
properties as linearly polarised electromagnetic waves in vacuum. 
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4.2 Alfven waves 

From © one obtains the equations of an Alfven simple wave: 

dB 2 dB 3 B 1 dE 1 dE 2 dE 3 



(47) 



E 2 -^ a B 3 E 3 +n a B 2 Va+Ci ^a(E 3 + fl a B 2 ) fJ, a (E 2 - (J, a B 3 ) ' 

where /i a is either /z„ or fi~ . First, we show that 

B 2 - E 2 = const (48) 
and 

Ha = const. (49) 

Since both these properties are Lorentz invariant it is sufficient to show that they hold in one particular frame. In this case, 
the most convenient frame is the "fluid" frame, where E — and, hence, dE 2 = 0. Moreover, ( f47j ) ensures 

d(B 2 ) = d(B 3 E 2 - B 2 E 3 ) = 

and the results (^j49) are manifest. 

Next, equations (17]) can easily be integrated to obtain 

rja = E 3 + fi a B 2 = const, (50) 

Ca = E 2 - (j, a B 3 = const, (51) 

( a E 2 + i] a E 3 = const, (52) 

E i B i = const, (53) 

Since the normal component of electric field varies the Alfven waves of FFDE are not quite transverse. This is similar to 
what has been found for Alfven waves in relativistic MHD (Komissarov 1997). These waves can still be described as linearly 
polarised. If we introduce the vector 

then 

dEtLt a and dB t \\t a . (54) 

From this equation and ( |46| ) one can see that the polarisations of Alfven and electromagnetic waves propagating in the same 
direction are not, in general, mutually orthogonal. 

If Bi = then from ( p6| ) it follows that Alfven waves propagate with the drift velocity 

lM a = (Ex B)i/B 2 , 

and, therefore, must be somewhat similar to the contact wave of relativistic MHD. Their properties are most transparent in 
the wave frame, where /i = 0. From (^6|) one can see that in such frame 

Et = 0. (55) 
Using this result and ( pS?] ) one finds that 

(-Si, B 2 , 0, 0, 0)* and (Si, 0, S 2 , 0, 0)' (56) 

is a basis of the Alfven eigenspace. The first vector describes the familiar rotation of St. The second allows variation of the 
magnitude of S t , but requires 

B 2 — E\ = const. 

Equations (^,^) show that all waves of FFDE are linearly degenerate, that is 
r ■ Vjy^i = 0, 

and, therefore, this system does not allow the formation of shocks via wave steepening. However, discontinuities may be 
introduced via discontinuous initial data and /or boundary conditions . For linearly degenerate waves the shock equations 
coincide with jump equations of simple waves ( Boillat I982| Anile 1989 ). 

Finally, figure I shows the normal wave speed diagram in the "fluid" frame, where E = 0. Since in this frame 

/i a = ±cos#, 

where 8 is the angle between S and the wave vector, all curves on this diagram are circles. 
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Figure 1. The normal speed diagram in the "fluid" frame (E = 0). 
5 ID UPWIND NUMERICAL SCHEME 

Nowadays, the construction of upwind numerical schemes for hyperbolic conservation laws is a well established area of numer- 
ical analysis. Here we describe a simple one-dimensional scheme which is very similar to the one constructed by Komissarov 



( 1999) for RMHD. The space-time is considered to be flat. 



5.1 General Structure 

The one-dimensional conservation laws of FFDE can be written in the form 

!n + d^ = Q > (57) 

where the vector of conserved variables is 

Q = {Si, 5*2, S3, B 2 , B 3 ) , 
where S is the Poynting flux vector, and the vector of hyperbolic fluxes is 

P = (Tn , T12, T13, —E3, E2) . 

Notice that we use all three space components of the energy-momentum equation. In addition, we introduce the auxiliary 
vector of "primitive" variables 

U — (B2, B 3 , Ei, E2, E3) . 
These primitive variables can be converted into the conservative ones and vice versa via 

S = E x B, (58) 

E=-^SxB. (59) 

The second equation automatically ensures the degeneracy condition (p^). Unfortunately, this does not mean that the second 
condition, (|l8|), will also be satisfied and extra care would have to be taken in cases where B 2 — E 2 is close to zero. At the 
moment we do not have the means to ensure (|l7|). 

Let us define a regular Cartesian grid such that the ith cell is centered at Xi = ih and occupies the region between 
x i-i/2 = (i — l/2)h and x i+1 / 2 — (i + l/2)h, where h is the mesh spacing. Now suppose that we know the solution at 
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t — t n and we want to calculate it at a later time t n +\. We integrate equations (^) over the ith cell and from t = t„ to 
t = t n+ i = t n + At to get 

Qi,n+i - + X (P *+*."+* - Pi-^n+i) = 0- (60) 
Here 



x ij- 1 



is the mean value of Q in the ith cell at time t n and 

tn + l 

is the time average of the flux at the cell interface with x = Xi + ^h. 

In a first order Godunov-type scheme (Godunov 1959), it is assumed that at the beginning of each time step the solution 
is uniform within each cell. This implies initial discontinuities at the cell interfaces and the time averaged fluxes can therefore 
be found by solving the corresponding Riemann problems. In fact, the solution, U* +1 j 2 , of the problem at the interface with 
x = 2^+1/2, does not depend on time but only on the initial left, Ui, n , and right, Ui+i t „, states. Then the first order fluxes 
can be computed from 

P i+i , n+i = P(u* +i ). 

In order to achieve second order accuracy, we use the first order scheme to obtain the solution, Q i n+ i and hence U i n+ i 
at the half time, t n+ i = t n + At/2. We then use this to compute average gradients of primitive variables in each cell as follows 



= "T av f L 
i h V 



„ _At7 ijn+ i,AC7 i+ljn+ , 

t,n-\- T£ 



(f) 

where 

AC/ iin+ i = (tf i>n+ i - U^_ hn 

and av(a, b) is a non-linear averaging function whose purpose is to reduce the scheme to first order in space in the neighborhood 



of discontinuities. Here we will adopt the same averaging function as in (Falle 1991) 
(a 2 b + ab 2 ) 



(a 2 + b 2 ) 



if ab > 0, and a z + b l 0, 



av(a, b) — 

These gradients can now be used to set up the left and right states for the second order Riemann problems 



if ab < or a 2 + b 2 = 0, 



The solution, U* +1 / 2 = U(U l i+1 / 2 , ^[+1/2)1 to this Riemann problem allows us to compute the second order fluxes 

P i+hn+i = P(u* +i ), 

which are used to advance the solution through the full time step from t = t n to t = t„+i according to (B3). 



5.2 Riemann Solvers 

By a Riemann problem we mean the initial value problem for ( |57| ) with the initial conditions 

For a general hyperbolic system this problem has a self-similar solution which involves only shocks and centered rarefactions 
(e.g. Landau & Lifshitz 1959, Jeffrey & Taniuti 1964). Since all waves of FFDE are linearly degenerate, centered simple waves 
do not exist and, thus, the solution involves only discontinuities. 
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5.2.1 Linear Riemann Solver 

Since all waves are linearly degenerate we have 



U r = Ul + C-iTi, (62) 

i=l,4 

where r*i is the eigenvector of the upstream state of ith wave. For fast waves these states are already known, they are the left 
and the right states of the Riemann problem. Upstream states of Alfven waves are different and depend on the amplitudes 
of fast waves. In our linear Riemann solver we ignore this difference. Then (^) becomes a linear system of five equations for 
four unknowns, Li. To handle it, we ignore one of the equations. The corresponding four components of the solution at x 1 — 
are then found via 

The remaining component is determined using the degeneracy condition (|l7j). 
5.2.2 Exact Riemann Solver 

In most Riemann problems originated from numerical simulations the difference between left and right states is rather small 
and linear Riemann solvers are sufficiently accurate. However, for handling strong discontinuities one might require an exact 
solution of the Riemann problem. Thanks to the simplicity of FFDE, an exact Riemann solver can be constructed rather 
easily. 

First we consider the case where none of the states of the Riemann problem has multiple eigenvalues. From the results 
of Sec. 4 it follows that only fast waves change the value of p = B 2 — E 2 and only Alfven waves change the value of E\. This 
implies that the Riemann problem reduces to determining the values of p and E\ in the region bounded by Alfven waves, p 
and Ei. Given the left state, the values of E\ and p in the right state, and p and E\ the corresponding right state V r (p,Ei) 
can be found by successive application of the jump equations from left to right: 

m = [Ss] = [Si] = m = m_ = m m) 

-m o n f Cf nfvt -2/^(^ + 0)' 

for fast waves, and 

[B 2 ] _ [B 3 ] _ gi[gi] = [E2I = _[Esl = [p] 

(a Va Va + (Z ^a^a ^aC,a 

for Alfven waves. This gives us the system of nonlinear equations 

Vr(p,E 1 ) = U r . (66) 



(65) 



which can then be solved iteratively. Once p and E\ are found the resolved state is determined via ( p4||65| ) . 
In degenerate cases the problem is even easier. 

(i) If [it = fi£ in the right state of the Riemann problem, then the right fast and Alfven waves merge into a single degenerate 
wave across which both p and E\ are constant Thus we immediately obtain p = p r , E\ — E\ T and then determine the resolved 
state in the same way as in the non-degenerate case. Similarly, if [ij = [i~ in the left state then p = pi, E\ = Ei t . 

(ii) It is easy to verify that for any of the degenerate waves, the corresponding degeneracy condition holds on both sides of 
the wave. Thus, if fit equals to [it in the right state and fij equals to [la in the left state then in the resolved state conditions 
( [32]) and (^) are satisfied simultaneously. This immediately gives E = and B± = in the resolved state. 

(iii) If Bi — then the Alfven waves merge into one contact-like wave across which both p and [i a are constant. This leads 
to the following two equations: 

pt(£ fl )= P * r (£ fr ), (67) 

=/C (£/,), (68) 

where * indicates states adjacent to the contact wave. Simple manipulations reduce this system to a quadratic equation for 
£ e , (or Ce r ). It appears that only one of its solutions satisfies the condition (|l8[). 

5.3 Test calculations 

We have tested the scheme against exact analytic solutions for plane electromagnetic and Alfven waves including all degenerate 
cases. The exact solutions can easily be constructed using the results of Sec. 3 and 4. The profiles of these waves should remain 
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unchanged as they propagate. Figures |2|-|3| illustrate the results for some of the test calculations. In these figures the initial 
and the final exact solutions are shown by continuous lines and markers show the final numerical solution. 

Fast wave: Figure 2a shows the results for a fast wave propagating in the positive direction. The initial solution is 

Si = 1.0, B 3 = E 2 = 0.0, 

!1.0 if x 1 < -0.1 

1.0 + | (a; 1 +0.1) if - O.K x 1 < 0.1 
1.3 if x 1 > 0.1 

Alfven wave: Figure 2b shows the results for a non-degenerate Alfven wave propagating in the negative direction with the 
wave speed /i = —0.5. The initial solution is 

B[ = B' 2 = 1.0, E' 2 = E' 3 = 0.0, 

1.0 if x 1 < -0.1 

B' z ={ 1.0 + § (i 1 +0.1) if - O.K x 1 < 0.1 , 
1.3 if x 1 > 0.1 

where B' and E' are measured in the wave frame which is moving relative to the grid with speed /i = —0.5. 

Degenerate Alfven wave: Figure (^) shows the results for a rotational degenerate Alfven wave propagating to the right with 
H — 0.5. The initial solution is 

E' = 0, B[ = 0, 

B 2 = 2 cos (f), B'j, = 2 sin <f>, 

where 

{0.0 if x 1 < -0.1 

^(x 1 +0.1) if — 0.1 < a; 1 < 0.1 . 
f if x 1 > 0.1 

Three waves problem: In this last test the initial discontinuity at x = splits into two fast discontinuities and a stationary 
Alfven discontinuity. The initial solution is 

B = (1.0, 1.5, 3.5) E = (—1.0, —0.5, 0.5) if x 1 < 0, 
B = (1.0, 2.0, 2.(3)) E = (-1.5, 1.(3), -0.5) if x 1 > 0. 

As one can see, numerical diffusion smoothes out the sharp wave fronts of the exact solutions but otherwise the agreement 
is very good indeed. 

Not all Riemann problems with both the left and the right states satisfying the degeneracy conditions ( fnjjli| ) have 
solutions within the frame of degenerate electrodynamics. For example, the following problem does not seem have one 

B = (1.0, 1.0, 1.0) E = (0.0,0.5,-0.5) if x 1 < 0, 
B = (1.0, —1.0, —1.0) = (0.0, 0.5, -0.5) if x 1 > 0, 

as iterations of our exact Riemann solver do not converge. To make sure that this is not due to a fault in the solver we 
have considered a related problem where this initial discontinuity is substituted with a linear transition layer. Numerical 
simulations (see figure ^) show that within this layer B 2 — E 2 decreases in time toward zero leading to the breakdown of the 
condition 

B 2 -E 2 > 0. 

As B 2 — E 2 tends to zero the drift velocity of plasma 
v d = E x B/B 2 

tends to the speed of light suggesting that one may no longer neglect the inertia of plasma particles. 



6 CONCLUSIONS 

We have shown that the evolution equations of FFDE can be written in the form of conservation laws. One group of equations 
arise from the energy-momentum conservation. There are only two linear independent dynamic equations (components of one 
tensor equation) in this group. This is in contrast to MHD, where all four components of the energy-momentum equation are 
linearly independent, and vacuum electrodynamics, where in general F ^ does not have zero eigenvectors and the energy- 
momentum equation is reduced to 
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Figure 3. The propagation of a rotational degenerate Alfven wave. The initial and the final exact solutions are shown by continuous 
lines and markers show the final numerical solution. 



V„F^ = 0, 

which has three independent dynamic components. The one-dimensional system of FFDE with plane symmetry, in flat space- 
time, includes four linearly independent dynamic equations for four independent components of the electromagnetic field 
tensor (other equations simply require B\ =const.) This system is hyperbolic, though not strictly hyperbolic, with a pair of 
fast waves and a pair of Alfven waves. All characteristic modes are linearly degenerate, which means that this system does 
not allow the formation of shocks via steepening of continuous waves. 
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Figure 4. In this test problem the initial discontinuity at x 1 = splits into two fast waves and one stationary Alfven wave. The exact 
solution at t = 0.75 is shown by continuous lines and markers show the corresponding numerical solution. 
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Figure 5. Breakdown of the approximation of degenerate electrodynamics as B 2 — E 2 — > 0. See the last problem described in Sec. 5. 4 



This formulation is particularly useful when it comes to developing numerical schemes for time-dependent FFDE as 
there has been a great deal of work on numerical methods for hyperbolic conservation laws. Here we have presented a 
simple one-dimensional scheme based mainly on our linear Riemann solver. Its generalization to multi-dimensions is relatively 
straightforward. We have already constructed a 2D numerical scheme adapted to Kerr space-time which will be described 
elsewhere. Such schemes can be useful tools in studying electromagnetically driven jets and winds from black holes and neutron 
stars and, perhaps, other related phenomena of relativistic astrophysics. 
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